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ABSTRACT 


In  this  thesis,  we  study  the  new  imaging  theory  developed  by  Cheney  and  Borden 
[1]  that  ineorporates  target  motion  during  data  collection  for  the  imaging  process.  The 
subject  of  radar  imaging  from  scattered  waves  is  explored  and  incorporated  into  the  new 
imaging  approach.  A  simulation  model  using  MATLAB  is  developed  to  simulate  the 
imaging  algorithm  and  also  to  validate  the  performance.  It  is  shown  that  the  new  imaging 
scheme  is  well  behaved  and  is  linear  shift  invariant  when  the  data  are  ideal.  It  is  also 
shown  that  the  geometry  of  the  transmitters  and  receivers  affects  the  behavior  of  the 
imaging  system. 
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I.  INTRODUCTION 


A,  RADAR  IMAGING 

Radar  is  a  system  that  uses  eleetromagnetie  waves  for  deteeting,  loeating,  and 
identifying  refleeting  objeets  over  long  distanees  in  both  adverse  and  good  weather 
eonditions.  Unlike  optieal  systems,  radar  systems  are  generally  not  affeeted  by 
atmospherie  attenuation  beeause  of  the  wavelengths  used.  However,  in  imaging,  the 
resolution  of  the  image  is  dependent  on  the  signal  wavelength.  Although  radar-based 
imaging  systems  have  eoarser  resolution  than  other  types  of  imaging  systems,  many 
objeets  of  interest  are  still  within  the  ranges  of  the  signal  wavelength  employed  by  radar 
systems. 

Many  radar  imaging  teehniques  have  been  developed  sinee  the  1950s  to  improve 
the  resolution  of  radar-based  imaging  systems.  Examples  of  sueh  teehniques  are 
Synthetie  Aperture  Radar  (SAR)  and  the  Inverse  Synthetie  Aperture  Radar  (ISAR), 
whieh  have  the  ability  to  produee  imagery  with  high  resolution. 

B,  IMAGING  TECHNIQUES  TO  CREATE  ARTIFACT  FREE  IMAGERY 

The  main  goal  in  radar-based  imaging  teehniques  is  to  ereate  an  artifaet  free 
image.  Model  assumptions  are  used  for  traditional  imaging  teehniques;  however,  these 
assumptions  indireetly  eause  the  image  artifaets  as  diseussed  in  the  literature.  The  start- 
stop  approximation  is  eommon  to  all  existing  imaging  teehniques  and  this  is  eentral  to  the 
thesis  and  holds  a  fundamental  differenee  between  the  new  imaging  seheme  developed  in 
[1]  and  eurrent  teehniques. 

1.  The  Born  Approximation 

In  the  standard  radar  seattering  model,  an  objeet  is  assumed  to  be  eomposed  of  a 
eolleetion  of  simple  point  seatterers.  When  eonsidering  multiple  seattering  (in  whieh  the 
seatterers  are  allowed  to  interaet)  analysis  beeomes  too  unwieldy.  Henee  the  “weak 
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scatterer”  or  Bom  approximation  is  used  instead.  This  approximation  treats  the  seatterers 
as  non-interaeting  and  results  in  a  more  manageable  linear  problem.  However,  with  the 
Born  approximation,  image  artifaets  may  still  exist,  as  this  approximation  diseards  the 
higher  order  multiple-seattering  terms. 

2.  Bandwidth  Limited  Radar  Systems 

Although,  ideally,  an  infinite  bandwidth  is  desired,  all  praetieal  systems  are 
bandwidth-limited,  and  this  faet  affeets  image  resolution.  In  addition,  the  signal 
information  measured  from  praetieal  systems  is  not  noise-free,  as  it  is  always  eormpted 
by  unwanted  signals,  and  this  also  eauses  unwanted  image  artifaets. 

3,  The  “Start-Stop”  Approximation 

Most  modern  radar  systems  use  a  train  of  pulses  together  with  eoherent 
integration  for  the  deteetion  of  a  target  and  for  an  estimation  of  the  target’s  veloeity. 
These  types  of  waveforms  typieally  allow  the  use  of  the  start-stop  approximation,  whieh 
assumes  the  target  is  stationary  during  the  measurement  proeess.  This  is  a  reasonable 
approximation,  sinee  the  speed  of  the  target  is  relatively  slow  when  eompared  to  the 
speed  of  the  eleetromagnetie  waves,  whieh  travel  at  the  speed  of  light.  By  this 
approximation,  veloeity  estimates  are  based  on  range  ehange-rate.  Target  position  and 
elosing  speed  ean  be  estimated  from  the  measured  range  and  range-rate  obtained  during 
the  measurement  proeess. 

In  order  to  aehieve  better  signal-to-noise  ratio  (SNR)  in  the  measurements, 
eoherent  integration  is  required.  However  this  ean  require  signifieant  number  of  pulses  in 
the  pulse  train.  When  the  pulse  train  gets  too  long,  the  stationary  target  approximation 
may  not  be  valid.  The  target  may  move  eonsiderably  within  the  pulse  interval  of  the  pulse 
train,  whieh  may  result  in  a  blurred  range. 
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C.  MOTIVATION:  IMAGING  MOVING  TARGETS 

Many  imaging  techniques  have  been  developed  to  image  moving  targets.  For 
example  the  Space-Time  Adaptive  Processing  (STAP)  (which  is  a  signal  processing 
technique)  uses  multiple-element  antenna  arrays  coupled  with  real-aperture  imaging 
techniques  to  produce  ground  moving  target  indicator  (GMTI)  images  [2],  The  technique 
described  in  [3]  uses  SAR  (designed  to  image  stationary  scenes)  together  with  GMTI 
processing  for  detecting  slow-moving  surface  targets  that  exhibit  start-stop  like 
maneuvers.  Another  method  described  in  [4]  discusses  an  implementation  based  on 
forward  looking  array  radar. 

Velocity  Synthetic  Aperture  Radar  (VSAR)  is  another  approach  and  is  described 
in  [5].  VSAR  is  a  multi-element  SAR  system  that  uses  conventional  processing  to  form 
an  image  for  each  element  (or  sub-array).  The  image  phases  are  preserved  and  compared 
to  estimate  the  target  velocity.  However,  this  method  assumes  that  the  scatterers  remain 
in  a  resolution  cell  throughout  the  measurement  period. 

There  is  a  common  assumption  made  by  these  imaging  techniques:  namely  that 
the  target  is  instantaneously  stationary  while  it  interacts  with  the  wave.  Therefore,  there  is 
generally  some  need  for  an  imaging  technique  that  can  accommodate  target  motion  in  a 
more  complete  way.  This  is  the  main  focus  of  this  thesis. 

D,  OBJECTIVE 

The  objective  of  this  thesis  is  to  study  a  linearized  imaging  theory,  developed  by 
Professors  Cheney  and  Borden  [1],  for  imaging  a  scene  with  moving  targets.  The  physics 
behind  it  and  the  approach  to  address  image  artifacts  associated  with  targets  moving  in  an 
unknown  fashion  will  be  discussed.  There  will  be  analysis  of  a  new  imaging  code  based 
on  a  moving  point  scatterer  via  MATLAB  simulation.  The  behavior  of  the  imaging  model 
will  also  be  examined  under  different  data  collection  geometries. 
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II.  IMAGING  THEORY 


In  this  chapter,  we  discuss  how  a  radar  system  obtains  range  and  Doppler 
information  of  a  target  and  how  radar  images  are  formed.  We  also  discuss  the  problems 
associated  with  radar  imaging  and  its  challenges. 

A,  RADAR  SYSTEMS 

Radar  systems  use  eleetromagnetie  waves  to  deteet  the  presence  of  a  target  in  an 
area  of  interest.  These  systems  send  out  a  signal,  whieh  could  be  a  pulse  or  a  series  of 
pulses  (known  as  pulse  train),  and  this  signal  interacts  with  a  target  and  is  refleeted  back 
to  the  radar.  From  the  reflected  signal,  the  two  measurements  that  ean  be  obtained  are  the 
round  trip  time  delay  of  the  transmitted  signal  and  the  target  radial  velocity.  The  range  of 
the  target  is  determine  by  measuring  the  round  trip  time  delay,  r ,  where  Range, 
R  =  CT  !  2  .  The  radical  velocity  can  be  obtained  from  the  phase  change  of  the  reflected 
signal  eompared  to  the  transmitted  signal. 

The  down  range  resolution  of  the  radar  system  is  determined  by  the  pulse 
bandwidth.  The  transmitting  frequeneies  are  usually  20  GHz  and  below,  as  the 
atmospherie  attenuation  of  waves  with  these  frequencies  is  negligible.  Higher  frequeneies 
are  subjeet  to  higher  atmospherie  attenuation.  However,  there  are  gaps  at  some  ranges 
among  the  higher  frequeneies  that  have  lesser  attenuation  and  eould  be  exploited  for 
various  requirements.  Systems  with  large  bandwidth  are  eapable  of  produeing  a  finer 
down  range  resolution,  and  such  systems  are  typically  sought  for  imaging  or 
identification  purposes. 

B,  SCATTERING  OF  ELECTROMAGNETIC  WAVES 

Radar  information  is  based  on  the  reflected  eleetromagnetie  waves. 
Understanding  the  behavior  and  properties  of  such  waves  provides  important  information 
about  the  target. 
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From  the  discussion  in  [6],  a  one-dimensional  scattering  model  for  a  moving 
pec'  plate  in  free  space  is  commonly  used  as  an  illustration.  Assuming  the  waveform 
generator  produces  a  time  varying  waveform  s{t)  and  is  mixed  with  the  carrier  wave  of 

frequency  to  produce  / {t)  =  s{t)cos{wQt) ,  the  transmitted  field  is, 


E,^Sx,t)  =  Erf{t-xlc)  (2.1) 

With  the  transmitted  wave  propagating  along  the  x-axis,  we  expect  the  scatter  field  to  be 
propagating  in  the  opposite  direction,  which  is  along  the  -x-axis  and  can  be  expressed 
as, 

E,,^,X^,t)  =  Erg{t  +  xlc)  (2.2) 

Further  assuming  the  position  of  the  plate  is  x  =  i?(t),  the  boundary  condition  on  the  PEC 
for  the  scattered  field  is, 

+  R{t)  /  c)  =  -E,„J{t  -R{t)/  c)  (2.3) 

For  simplicity,  assuming  that  the  plate  undergoes  linear  motion  so 
thati?(t)  =  x  +  vt  ^^^(^tboundaiy)  ’  whcrc  V  is  the  rate  the  range  changes  and  is  also  known  as 

the  range-rate.  For  short  pulses,  R{t)  =  x  +  vt  can  be  considered  to  be  the  first  two  terms 
of  Taylor  expansion  (valid  for  small  t).  Substituting,  u  =  t  +  R{t)  /  c ,  wq  can  solve  for  t 
in  terms  of  u  (at  boundary  condition): 


t  = 


u-R!  c 
\  +  v  I  c 


(2.4) 


and  so  the  scattered  field  can  be  expressed  as: 


Escatt  g{t  +  R{t)l c)  =  -Einc  f\a{t  +  xl c-R! c)-R I c\  (2.5) 

where  the  Doppler  scale  factor  a  is: 

^  A  perfect  electrical  conductor  (PEC)  allows  the  charges  to  move  freely  and  instantaneously  in 
response  to  a  field  and  the  fields  inside  a  PEC  are  zero.  The  boundary  conditions  are  that  the  tangential 
components  of  the  electric  field  musf  be  zero  and  the  tangential  components  of  the  magnetic  field  are 
relafed  to  currents  flowing  on  the  surface  of  the  PEC. 
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(2.6) 


1-v/c 

a  = - 

1  +  v/c 

Since  the  transmitted  waveform  is  f(t)  =  s(t)  cosico^t)  ,  the  scattered  signal  obtained  at  the 
antenna  (x  =  0)  is: 

Prec  (0  =  /[<^(^  -  R/  c)-  R/  c]  =  s[a{t  -R/  c)-R/  c]  cos[Wo  {a{t  -  R  /  c)- R  /  c)]  (2.7) 

for  a  scatterer  that  is  slowly  varying  so  that  v  /  c  is  small.  We  can  expand  the 
denominator  of  Equation  (2.6)  in  a  geometric  series  to  obtain 

«  =  (l-v/c)(l-v/c  +  0[(v/c)']  =  l-2v/c  +  0[(v/c)"]  (2.8) 

approximating  a  »  1  -  2v  /  c  ,  we  can  approximate  the  argument  a{t  -  R!  c)-  R!  c  ~ 
t  -2RI  c  to  obtain 

Precis)  ~  s{t-2R/ c)cos(oq  [(1-2v/ c){t-R/ c)-R/ c])  (2.9) 

It  can  be  seen  that  the  carrier  frequency  has  been  shifted  by  an  amount, 

2v 

C0j^= - cOq  (2.10) 

c 

This  is  known  as  the  Doppler  shift.  If  the  target  moves  towards  the  antenna  (v  <  0),  then 
the  Doppler  shift  is  positive.  If  the  target  moves  away  from  the  antenna  (v  >  0),  then  the 
Doppler  shift  is  negative.  Although  radar  targets  are  not  usually  made  of  PEC  plates,  the 
received  signal  is  nevertheless  well  modeled  as  a  time  delayed  and  Doppler-shifted 
version  of  the  transmitted  signal.  This  model  is  valid  for  most  targets. 

C  CORRELATION  RECEPTION 

Erom  the  radar  range  equation  [7],  it  is  easily  seen  that  the  received  signal  is 
reduced  in  energy  by  a  factor  of  R  of  the  transmitted  signal.  Earger  detection  ranges  can 
be  aehieved  by  inereasing  the  transmitted  power  of  the  signal;  however  the  transmitter 
power  is  usually  limited  by  operational  requirements.  It  is  common  for  the  receiver  of  the 
radar  to  be  expected  to  detect  a  signal  power  as  low  as  10'  Watts. 
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Although  we  can  easily  increase  the  signal  to  noise  ratio  (SNR)  by  doing  coherent 
pulse  integration,  there  are  limitations.  For  example,  the  increase  of  pulse  repetition 
frequency  will  also  decrease  the  maximum  unambiguous  range.  The  target  usually  moves 
during  data  collection,  which  results  in  the  phase  of  the  scattered  field  being  altered. 
Pulse  integration  is  also  limited  by  the  ability  of  the  local  oscillator  to  remain  coherent 
over  the  interval  in  which  the  pulses  are  transmitted. 

Radar  obtains  information  about  a  target  by  comparing  the  received  signal  with 
the  transmitted  signal.  In  such  processing,  the  objective  is  to  separate  the  scattered  field 
from  the  random  noise  contamination  in  the  received  signal  (t) . 

Srecit)  =  S,,,nit)  +  n{t)  (2.11) 

Maximum  likelihood  processing  selects  the  best  estimate  for  the  scattered  field,  which 
maximizes  the  probability  of  scattered  signal  for  a  given  over  all  time. 

From  [6],  it  is  shown  that  such  processing,  averaged  over  all  time  reduces  to  the  problem 
of  finding  a  function  that  minimizes: 

00  00  00 

i  l\„(Or*'+  i  2 Re  I  (<>;„('')*'  (2.12) 

—00  —00  —00 

Searching  for  the  best  minimizing  function  is  made  more  efficient  by  restricting 
the  search  to  a  few  parameters.  In  this  case,  the  natural  model  is  based  on  the  scattering 
interaction  between  the  interrogating  field  and  the  target.  If  denotes  an  incident 

pulse  transmitted  by  the  radar,  then  the  linear  radar  (weak-scatterer)  scattering  model 
follows  by  superposition 

00 

=  (2.13) 

-00 

where  p{o,t)  is  the  target  reflectivity  density  defined  such  that  p{u,T)dTdu  is 
proportional  to  the  field  reflected  from  the  target  at  range  between  cr  /  2  and 
c(r  +  (ir)/ 2  with  Doppler  shift  between  u  and  u  +  du. 
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From  Equation  (2.13)  ,  the  target  refleetivity  seattered  field  is  dependent  on  the 
two  parameters  r  and  o .  The  parameters  for  the  simplest  seareh  spaee  would  be  the  time 
shift  T  and  the  Doppler  shift  o .  The  eorrelation  reeeiver  tries  to  find  r  and  o  whieh 
maximize  the  real  part  of: 

00 

ri{o,T)  =  J  (2.14) 


The  real  part  of  the  eorrelation  integral  output  will  be  maximum  when  the  reeeived  signal 
“looks-like”  a  time-delayed  and  frequeney-shifted  version  of  the  transmitted  signal.  Other 
signals,  like  noise  from  the  reeeiver  eireuit  and  the  environment,  will  be  suppressed.  As 
sueh,  the  eorrelation  reeeiver  appears  to  be  operating  like  a  type  of  “matehed  filter”. 

By  substituting  Equation  (2.13)  into  Equation  (2.14),  the  output  of  the  eorrelation 
reeeiver  results  in 


00 

ri{v,T)  =  \\  p(u',t')z(o-u',t 

-00 


T')e^  dr'do' 


(2.15) 


where 


Zio,T)=  j 


(2.16) 


The  eorrelation  noise  term  is  eonsidered  to  be  very  small  as  eompared  to  the  signal  of 
interest,  and  henee  it  is  not  explieitly  displayed  in  the  Equation.  The  standard  radar  data 
model  is  shown  in  Equation  (2.15),  and  it  expresses  the  output  of  the  eorrelation  reeeivers 
as  a  eonvolution  of  p  and  z  ■  The  function  j(t>',r)  as  defined  in  Equation  (2.16)  is  also 
known  as  the  ambiguity  function. 


D,  ONE-DIMENSIONAL  (HIGH  RANGE  RESOLUTION)  IMAGING 

In  general,  radar  pulses  are  often  of  short  duration  and  target  velocity  is  almost 
stationary  when  compared  to  the  speed  of  light.  These  short  pulse  signals  can  be  used  to 
form  one-dimensional  high  range  resolution  (HRR)  “images”  of  the  target.  This  method 
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works  when  the  transmitted  pulse’s  instantaneous  range  resolution  is  smaller  than  the  size 
of  the  target.  When  the  seattering  of  the  sub-elements  of  the  target  are  non-interaeting 
and  point-like,  the  scattered  pulse  will  be  a  sum  of  the  damped  and  blurred  images  of  the 
incident  pulse,  which  are  shifted  by  time-delays  that  are  proportional  to  the  sub-element’s 
range  [8].  Figure  1  shows  a  one-dimensional  image  created  by  an  HRR  radar  system.  It  is 
clearly  shown  that  the  peak  profiles  are  associated  with  isolated  target  elements. 


y  (motoi*s) 

0  10  20  M)  -40 


Figure  1 .  Example  of  a  range  profile  from  a  B-727  jetliner.  The  top  view  (with 
orientation  at  time  of  measurement)  is  displayed  beneath  (From:  [8]) 


The  received  signal  of  non-interacting,  stationary  targets  can  be  expressed  as 

00 

Srec(0=  J  PiT')Si„,{t-T')dT'+n{t)  (2.17) 

-00 

where  n{t)  is  the  random  noise  term.  The  correlation  integral  can  be  expressed  as 
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vi^)=  I  Srecit'yincit'-^)dt'  =  J  4c p{t ')s '- T ')dT ' dt '  +  Hoisc  term 

—00  —00  —00 

00 

/7(4  =  jJ  sl^(t'-'i^)Si„^(t'-T')dt' p(T')dT'  +  noise  term  (2.18) 

-00 

Substituting  t"  =  t'-  p  and  dropping  the  noise  term  gives 

00 

?7(4=  j  X{^-^')p{^')dr'  (2.19) 

-00 

where  x  is  the  autocorrelation  function  (in  this  case,  it  is  also  a  point  spread  function.) 

00  00 

X'(r)  =  J  5*(t"-r)5(t")t/t"  =  1  s* {t')s{t'+r)dt'  (2.20) 

-00  -00 

The  image  is  formed  by  the  convolution  of  p  and  x, 

However,  range  profiles  can  be  difficult  to  be  used  for  target  classification,  since 
all  scatterers  located  at  the  same  distance  from  the  radar  will  reflect  the  signal  back  with 
the  same  time-delay.  As  a  result,  the  range  profile  will  not  be  able  to  distinguish  the  cross 
range  structure. 

E,  TWO-DIMENSIONAL  IMAGING 

A  simple  concept  based  on  “triangulation”  can  be  used  to  determine  the  cross¬ 
range  target  structure  while  still  using  an  HRR  radar  system.  Figure  2.  illustrates  this 
approach,  which  relies  on  range  profiles  collected  from  different  target  orientations  and 
correlated  to  form  a  two-dimensional  image. 

Consider  a  set  of  three  point  targets  with  the  radar  located  at  the  same  distance 
from  targets  2  and  3  in  Figure  2.  Given  the  orientation  as  shown,  the  return  signal  will 
only  indicate  two  targets.  Hence,  ambiguity  exists  when  targets  are  the  same  distance 
away  from  the  radar. 
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Target  1 

o 


Bands  of  constant 
range  from  radar 


o 

Target  2 


Target  3 

o 

Single  i)ulse  ainbignity 


Figure  2.  Ambiguous  scenario  from  a  single  radar  pulse  (From;  [6]) 


The  return  signals  from  multiple  directions  can  be  used  to  determine  range  and 
cross-range  information.  Hence,  the  range  profiles  are  swept  in  the  cross  range  direction 
to  form  bands  representing  possible  locations  of  target  scattering  centers  in  space.  The 
bands  are  then  superimposed,  and  the  crossing  points  are  used  to  determine  the  scattering 
center  locations. 

Multiple  inilse  ‘■imaging” 

Target  1 

o 

Target  3 

t)  o 

Target  2 

I  A 

Figure  3.  Cross-range  information  obtained  from  range  profiles  (From;  [6]) 


There  are  2  different  schemes  used  to  collect  data  from  different  target  aspects.  In 
one,  the  radar  moves  while  the  target  remains  stationary.  This  data  collection  process 
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occurs  over  a  synthetic  aperture  and  henee  is  called  “synthetic  aperture  radar”  (SAR).  In 
the  other,  the  target  rotates  while  the  radar  is  stationary  and  staring  at  the  target.  This  data 
eolleetion  proeess  is  called  “inverse  synthetie  aperture  radar”  (ISAR). 

F.  RADAR  IMAGING— AN  INVERSE  PROBLEM 

Normal  radar  can  be  approximated  as  a  linear  measurement  system,  which  is  also 
shift  invariant.  Radar  imaging  would  not  be  possible  if  the  time-delayed  seattering  signal 
looked  different  for  different  time  delays.  Hence,  a  radar  imaging  system  can  be 
expressed  by  the  functional  operator  or  “kernel”  k  that  describes  how  the  measurement 
system  works 

m  =  Kf  (2.21) 

Where  m  denotes  the  measurement  at  position  x  and  /  denotes  the  “object”  function. 

The  direct  problem  in  radar  imaging  refers  to  predicting  the  data  measured  by  the 
radar  from  a  known  target.  The  inverse  problem  is  the  “reproduction”  of  the  original 
unknown  target  from  given  radar  measured  data. 

1.  Well-Posed  and  Ill-Posed  Problems 

Radar  systems,  whieh  are  typieally  band-limited,  do  not  allow  a  unique  solution  to 
the  inverse  problem.  This  is  due  to  the  faet  they  yield  finite  measurement  data,  which  are 
noisy  as  well. 

Using  diserete  measurement  systems  in  whieh  the  measurement  data  are  finite. 
Equation  2.16  can  be  written  in  matrix  form.  For  such  systems,  the  goal  is  to  find  a  good 
estimate  /  of  /  from  the  measurements  m  which  can  be  expressed  symbolically  as 

7  =  X-  m  (2.22) 

The  fact  that  the  “kernel”  k  is  not  always  a  square  matrix  poses  some  problems.  In 
general,  the  matrix  k  is  N  x  M,  and  when  M  >  N  there  are  more  unknowns  than  there 
are  linearly  independent  Equations  and  so  the  system  cannot  have  a  unique  solution. 
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As  the  object  space  (M-dimensional)  is  larger  than  the  measurement  space  (N- 
dimensional),  a  null  space  will  exist,  and  it  will  consists  of  all  vectors  of  the 
form  /-/c/.  The  null  space  is  determined  by  the  kernel /c ,  which  represents  the 
measurement  system.  As  such,  the  null  space  contains  all  things  that  cannot  be  measured. 
Null  space  accounts  for  image  artifacts  and  must  be  carefully  studied  when  interpreting 
reconstructed  images. 

In  order  to  solve  the  inverse  problem,  it  is  often  sufficient  to  search  for 
approximate  solutions  satisfying  additional  constraints  [9].  Since  the  set  of  approximate 
solutions  is  the  set  of  objects  with  images  closest  to  the  measured  one,  /  must  obey 

K^Kf  =  K^m  (2.23) 

This  matrix  Equation  is  also  known  as  the  “Normal  Equation”. 

2,  Data  Reconstruction — Regularization 

Equation  (2.23)  is  the  set  of  conditions  that  must  be  satisfied  by  the  least-squares 
solution  /  which  can  be  expressed  as 

f  =  {K^Ky^K^m  (2.24) 

Two  techniques  are  commonly  used,  namely  the  singular  value  decomposition  (SVD)  and 
Tikhonov  regularization.  Erom  linear  algebra  methods,  the  matrix  {k^kY^k^  can  be 
expressed  in  terms  of  its  SVD  [6]; 

{/c^K-y'/c^  =UDV^  (2.25) 

where  the  matrices  U  and  V  contain  the  eigenvectors  of  and  D  contains  the 

reciprocal  of  the  eigenvalues,  which  has  diagonal  elements  of  the  form  1  /  (where 

are  the  corresponding,  non  zero  eigenvalues).  Since  k  is  bounded,  the  eigenvalues  must 
form  a  sequence  that  gets  close  to  zero  as  the  number  of  measurements  (dimensions) 
increase,  so  their  reciprocals  become  arbitrarily  large.  Since  the  data  m  are  contaminated 

by  noise,  then  the  solution  will  multiply  the  noise  by  arbitrarily  large  values.  Hence,  the 

14 


simplest  approach  to  mitigate  the  noise  effect  is  to  truncate  the  SVD  representation  and 
replace  the  component  of  D  by  zero  if  the  corresponding  eigenvalue  falls  below  some 
threshold  value.  This  method  is  called  the  truncation  filter. 

Tikhonov  regularization  makes  an  ad-hoc  modification  to  the  “Normal  Equation” 
and  can  be  rewritten  as 


/  =  {K^K  +  aI) 


(2.26) 


Where  a  is  a  fixed  scalar.  When  a  =  0 ,  Equation  (2.26)  reverts  to  the  “Normal 
Equation”  and  so  it  can  be  expected  that  small  values  of  a  will  not  change  the  system 
description  significantly. 

After  introducing « ,  the  diagonal  elements  of  D  in  Equation  (2.25)  can  be 
expressed  as 


A.  = 


_  A 


X  +  a 


(2.27) 


when  X  »  ot ,  D..  «  1  /  y X.  (as  before).  However,  when  X.=  a,  D..  «  /  a  and  will 

not  become  arbitrarily  large  and  the  noise  issue  is  moderated. 

Eor  limiting  the  effect  of  noise,  both  Truncation  filters  and  Tikhonov 
regularization  are  powerful  and  useful  methods.  However,  none  offers  a  method  to 
choose  the  “best”  threshold  value  a  .  With  too  large  a  ,  the  ability  to  accurately  estimate 
the  object  will  be  affected,  and  too  small  a  will  not  control  the  noise  effectively. 
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III.  IMAGING  AND  THE  MOVING  TARGET  DATA  MODEL 


The  previous  chapter  discussed  different  types  of  imaging  theories  and  techniques 
that  are  used  for  radar  imaging.  This  chapter  describes  the  ideas  and  concepts  used  to 
develop  the  imaging  model  for  moving  targets. 

A.  LINEARIZED  DATA  MODEL 

The  scattering  model  used  to  describe  non  moving  targets  is  based  on  a  scattered 
field  using  Born  approximation  and  is 

scan  0  =  g(x x;  t t)p{x  (x t  ')d'x '  (3.1) 

where 


g(x\x;t',t) 


d{t-t'- 1  x'-x  \  Ic) 
An\x'-x\ 


(3.2) 


is  the  Green’s  function  and,  for  a  signal  5  ,(t)  transmitted  from  position  y  with  starting 

y 

time  -T  ,  it  can  be  shown  that  the  incident  field  at  x  is 

V  ’ 


Sy{t  +  Ty-\x-y\lc) 
A7:\x-y\ 


(3.3) 


The  scattered  field  WscatXy-'^d)  at  position  z  and  time  t  is  therefore 


¥scatXy^^d) 


z  — x' 

|/c)  Sy{t'+Ty- 

- 

x'-. 

y 

/c) 

4;r| 

z  — x' 

An 

X 

■’-Tl 

■p{x')d^x'  dt' 


(3.4) 


By  using  the  Born  approximation  and  the  time-domain  Green  function,  we  obtain  the 
scattered  field  for  non-moving  scatterers  as  shown  in  Equation  (3.4).  However,  we  are 
generally  interested  in  the  scattered  field  from  moving  scatterers  (targets),  hence  the 
scattering  model  shown  in  Equation  (3.4)  has  to  be  modified. 
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The  scattering  density  function  p{x)  used  in  Equation  (3.1)  is  only  space 
dependent,  since  the  target  is  not  moving.  It  is  like  a  scaling  function,  which  attenuates 
the  transmitted  signal.  For  a  moving  target,  the  scattering  density  function  is  dependent 
on  both  distance  and  time,  so  it  is  expressed  as  p(x,t) . 


So,  we  modify  the  scattering  model  by  letting  p^{x  -  vt)d^xd\  represent  the 

scatterers  in  the  volume  element  d^xd^v  of  phase  space  centered  at  position  x  and 
velocity  v.  Choose  T  to  be  the  time  of  the  transmitted  signal  such  that  the  scattering 

density  p^  =  p^{x) ,  which  accounts  for  all  targets  that  are  moving  with  velocity  v  at  time 
t=0  (i.e.,  the  time  of  the  transmitted  signal  is  activated  at  time  t  =  -T  ).  This  shows  that 
the  spatial  scatterer  density  centered  at  time  t  and  position  x  is 

=  (3-5) 

Hence,  the  scattered  field  is 


¥sca„iy,^d)  =  \\\ 


^(t-t’-|z-x’|/c)^^fy(C+fy-|x’-y|/c) 

4;r  I  z  -  X '  I 


p^{x'-vt')d^vd^x'  dt'  (3.6) 


4;r  I  x'-y  I 

By  implementing  the  change  of  variables  x'  ^  x  =  x'-vt' ,  the  modified  scattered  field  is 
■S(t-t'-\x  +  vt'-z\/c)  Sy{t'+Ty-\x  +  vt'-y\/c) 


¥scatXy^^d)  =  \\\- 


4;r  I  X  +  vt z 


4;r  I  X  +  vt y  I 


- p^{x)d  vd  xdt'  (3.7) 


Now,  the  scattered  field  of  moving  targets  is  incorporated  into  the  data  model.  The 
physical  interpretation  of  Equation  (3.7)  is  that  the  electromagnetic  wave  generated  from 
transmitter  at  position y  at  time  -T  interacts  with  the  target  at  time  t' .  During  the 

interval  [0,  t '] ,  the  target  has  moved  from  position  x  to  a  new  position  x  +  vf.  The 
scattered  field  with  strength  p^{x)  propagates  from  its  new  position  x  +  vt'  to  receiver  at 
location  z  and  at  time  t . 

Equation  (3.7)  can  be  simplified,  by  letting 
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R  =  R/R 


R^  =  x  +  vt  -  z,  i?=|/f|, 


and  substituting  into  the  scattered  field  Equation  (3.6)  we  obtain 


VscatXy^^^t)  =  \\\  - ,  „  - X - ,  „  - PX^)d  vd  xdt 


(3.8) 


The  dependence  of  R^^{t')  on  the  variable  t' makes  this  integration  complicated.  To 
further  simplify,  we  let  t'  =  ^{t)  ,  which  denotes  the  implicit  solution  of 

t-t'-  R^Xt')!  c  =  0,  i.e.,  t-JXt)-RX\Xt))lc  =  0 


The  integration  over  t '  then  yields 


/  3  ff^3.(4.v(0  +  7;,-^.,^(4.v(0)/c) 

¥scattiy^^X)  =  - ^ - = - = - - p  {x)d  vd  X 

JJ  {Anf  R^XKMR,^^{t,Xt)) 


(3.9) 


By  substituting  t  {t)  =  t-R  {t  (t)) / c ,  we  obtain 


¥scm  Ky,^X)=  — ^ - 3 - = - = - ^  Pv  {x)d  vd  X  (3.10) 


B,  FURTHER  SIMPLIFYING  APPROXIMATIONS 

Equation  (3.10)  is  applicable  to  most  general  scenarios.  It  is  applicable  to  fast 
moving  targets,  where  the  start-stop  approximation  is  not  valid.  It  can  also  be  used  for 
any  types  of  transmitted  waveforms. 

However,  retarded-time  problems  are  still  present.  These  problems  arise  due  to 
the  definition  of  the  retarded-time  ^,(0  .  The  retarded-time  can  be  obtained  by  solving 

4,v(0  =  ^-^x,z(4,v(0)/c  (3.11) 
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Hence,  a  more  meaningful  result  can  be  obtained  when  the  scatterers  are  moving  in  a 
known  fashion  so  that  R^^{t)  can  be  expanded  in  a  Taylor  series  around  t  =  0  and 

approximated  by  retaining  the  terms  linear  in  t.  We  refer  to  this  approximation  as  “slow- 
moving”  target. 


The  Slow-Mover  Approximation 


Assuming  that  (|  v  1 1)  and  (|  v  |  r  x  /  c)  are  much  less  than  \x-  z  \  and 
\x-  y\,  where  co  denotes  the  maximum  angular  frequency  of  the  transmitted 


signals ,,  we  have 


^xz(0  =1  Z-{X-Vt)  |=i?^^(0)  +  ifx,z(0)-Vt  +  ... 


(3.12) 


where  if^^(0)  =  jc-z,  =|  | ,  and  /fx,z(0)  = .  Substituting 

this  result  into  the  definition  of  retarded  time  yields, 


4,v  (0  ~  ^  -  ^x.z  (0)  +  (0)  •  (0 


/  C  ; 


t-K,Mic 

l-i-ifx,z(0)-v/  c 


(3.13) 


Inserting  this  approximation  into  the  result  for  and  simplifying  yields 

Wscatt(y’^’^)=  [[ - ^ - -p^{x)d  vd  X  (3.14) 

"  JJ  K,MK,yiO) 

where 


a  = 

x,v 


\-Rx,y{0)-V  /  C 
l  +  i?x,z(0)-v/  C 


(3.15) 


^  is  the  Doppler  scale  factor. 
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2. 


The  Slow-Mover  and  Narrow-Band  Approximation 


Radar  systems  are  typieally  narrow-band,  and  the  transmitted  signal  is  usually  of 
the  form 

Sy{t)  =  Sy{t)e^‘"’^‘  (3.16) 

where  Sy{t)  is  slowly  varying  when  eompared  to  e  ,  with  the  earrier  frequeney  of 

the  transmitter  y .  The  seeond  time  derivatives  of  (t)  are  dominated  by  the  e 

faetor,  and  by  implementing  the  slowly  varying  and  narrow-band  approximation,  we 
obtain, 

2  i(p.  „  -ico„a^  J 

CP  CD  e  •  e  ^  • 

°  (3.17) 

(^  +  ^3'  “  (0)  +  K,y  (0))  /  c)  P,  {x)d\d^x 

where 

ft,.  =  <’>]  [^,.,(0)  -  cT;  +  a„,7!„(0)]  /  c  (3.18) 

3,  The  Slow-Mover,  Narrow-Band  and  Far-Field  Approximation 

When  the  distanees  from  the  transmitter  to  target  and  target  to  the  reeeiver  are 
large  in  eomparison  with  the  imaging  seene  dimensions,  then  |  x  +  vt  |  and 

\x  +  vt?  xcD  I  c  ean  be  assumed  to  be  mueh  less  than  either  I  z  I  or  I  y  I .  Therefore,  the 

expansion 

A 

=1  Z-{x-Vt)  1=1  Z  I  -z-{x  +  vt)  + ... 

ean  be  applied  (and  similarly  for  Substituting  the  expansions  into  previous 

approximations  Equation  (3.17)  then  yields 
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(3.19) 


y-x)l  cjp^  {x)d^vd^x 

where 

A  A 

(P.,v=^y  \y\-yx-cT^  +  aX\z\-z-x)  Ic  (3.20) 

The  Doppler  seale  factor  can  also  be  approximated  as 

A 

1- V- v/ c  .  y 

— X - ~\-{y+z)-vlc  (3.21) 

l  +  z-v/c 

A  A 

because  |  v  |  /c  « I .  The  quantity  co^  x  is  the  Doppler  shift  (where /3^  =  -{y+  z)-v I c) 
By  inserting  T  =\y\lc  and  /  c  wq  obtain 

(P.,v=k^\z\-kyiy+z)-  x  +  {z-{z-x))vl c  (3.22) 

C.  IMAGING  VIA  A  FILTERED  ADJOINT 

In  developing  the  correlation  receiver,  it  is  observed  that  how  much  two  signal 
functions  (complex-valued)  look  alike  is  given  by  the  cross  correlation.  The  correlation 
integral  measures  how  much  the  received  signal  “looks  like”  the  time-delayed, 
frequency-shifted  version  of  the  transmitted  signal.  A  large  value  of  the  correlation 
integral  indicates  a  strong  resemblance,  while  a  small  value  indicates  a  weak 
resemblance. 

Using  the  same  notion,  the  usual  radar  model  is  developed  for  the  scattered  field 

of  a  point  target,  which  is  a  time-delayed  and  Doppler  shifted  version  of  the  incident 

field.  This  two-parameter  model  is  a  cross  correlation  between  the  received  signal  and  the 

transmitted  signal.  The  values  of  the  time-delayed  and  Doppler  shift  parameters,  which 

maximize  the  cross-correlation,  are  the  parameters  that  most  represent  the  target. 
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(4^)  I  2  Ilk  I 


e  e  ^ 


t  +  Ty-{\z\-z-x+\y\- 


The  generalized  model  for  depends  on  the  position  of  the  scatterer 

X  =  {x^,x^,x^)  and  its  velocity  v  =  {v^,v^,v^)  .  Using  y/ to  build  a  six-parameter 
scattering  model  based  on  the  arbitrary  position  p  =  {p^,p^,p^)  and  arbitrary  velocity 
u  =  {u^,u^,u^)  of  any  scatterer,  the  values  of  p  and  u  that  maximize  the  cross 

correlation  between  the  arbitrary  model  and  the  measured  data  will  represent  the  “best 
fit”  of  p  and  u  to  the  data.  These  values  of  p  and  u  will  localize  the  position  and 
velocity  of  any  unknown  scatterers  in  the  position-velocity  space  (phase  space). 

The  scattered  field  from  known  scatterers  p^{x)  with  positions  x  and  velocities  v 
is 


_^ycoyt-\Alc) 

(47rf\z\\y\ 


II  exp  I  -iky  (y+  z)  -  [x  -v{t  -  z  ■  (z  -  x)  /  c)] 


(3.23) 


xs. 


t-{\z\-{y+z)-x)l c  p^{x)dvdx 


To  simplify,  substitute  f  =  t-\z\/c , 


¥scatXy^^d')'^ 


2  -i(oJ 

-oiye 

(4;r)fyz||y| 


II exp  [-iky  {y+  z)-[x- v(t '+z-x/  c)] 


(3.24) 


X5, 


t'+{y+z)-x/ c  p^{x)d^vd^x 


By  ignoring  the  intensity  prefactor,  a  parametric  model  representing  the  field  from  an 
unknown  point  source  is 


y/iy, z,t')  =  exp  -ik  (y+  z)-[p -u{t'+ z- p  !  c)] 


X5. 


A  A 

f+  {y+z)-plc 
V  y 


(3.25) 


Based  on  the  cross-correlation  notion,  an  image  can  be  created  as  follows: 

/(;?,  m)  =  III  {y,  z,  t')y/*{y,z,t')dt' d’”  yd"  z  (3 .26) 
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where  m ,  n  depends  on  the  number  of  the  transmitter(s)/receiver(s),  and  the  integrals  are 
over  all  values  of  t' ,  y  and  z  for  which  data  are  measured.  Large  values  of  \I{p,u)\ 
indicate  a  strong  resemblance  while  low  values  of  | |  indicate  a  mismatch.  True 
cross  correlation  will  not  be  possible  if  the  data  measured  is  not  for  all  t' ,  y  and  z , 
therefore  the  best  approximation  to  cross  correlation  with  limited  data  will  result  in  image 
artifacts.  Due  to  these  artifacts,  a  fdter  function  Q{co,t',p,u,y,z)  is  included  to  suppress 
those  artifacts, 

I{p,  u)  =  jjj  {y,  z,  t  V*  {y,  z,  t  t\p,u,y,  z)dt '  d'^yd’^z  (3 .27) 

Hence,  the  image  is  formed  as. 


I{p,u)  =  -  jjj  Qio),  t\p,u,y,  z)e““^‘  e 


iky(y+^)- 


p-u{t'+z-p/c) 


xs. 


t'+{y+z)- p/ c  {y,z,t ')dt 'd"'yd"z 


(3.28) 


D,  IMAGE  ANALYSIS 


To  determine  the  effectiveness  of  this  imaging  model,  we  work  various  scenarios. 
By  inserting  the  ideal  measured  data  back  into  the  imaging  Equation  (3.28)  and 
substituting  {y,  z,  t ')  we  obtain. 


JJJ  JJ  (4;r)^  I  ^  II  T  I 


Pv(^) 


xexpa'A:  (y+z)-[/?-x-(M  -v)t'-u{z- p)/ c  +  v(z-x)/ c]>  (3.29) 


t'+  (y+z)- p/ c jSy  '+ (_y+ z) •  X / c jd  vd  xdt' d"" yd" z 


The  distances  |  z  |  and  |  y  \  are  amplitude  scaling  factors  and  not  important  for 
imaging  purposes.  Thus  we  can  choose 


Qico,t\p,u,y,z)  = 


(4;r)"  |z||y 


0).. 


J{p,u,y,z) 


(3.30) 
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where  J(p,u,y,z)  is  geometry  dependent  and  is  used  to  eompensate  for  the  Jacobian  that 
results, 


=  t'+{y+z)-p/c  Sy  t'+{y+z)-x/c 


J  ^ 


-ik  (y+z)iu-v)t' 


X  y{y+  z)  -  {p  -  X  - u{z-  p)  I  c  +  v{z-x)l  c] 

xj(p,u,y,z)p^{x)dt’d\d^xd"'yd"z 
This  result  can  be  written  as 

I{p,u)  =  jj  K{p,u,y,z)pXx)d\d^x 

where 


(3.31) 


(3.32) 


K{p,u,y,  z)  =  jj exp  j -ik^ {y+  z)  ■  [u{z-  p)  -  v(z- x)]  /  c 


xj'?v[  ^'+(T+^)-.p/cJ^3,^t'+(y+z)-x/c 


X  J{p,  u,  y,  z)d'’'yd"z 


(3.33) 


is  the  point-spread  function  describing  the  behavior  of  the  imaging  system.  By  applying 

A  A  A  A 

the  change  of  variables  t  =  t'+  —  {y+ z)-{p  +  x) / c  and  set  r  =  {y+  z)-{p-x) / c  we 
obtain. 


K{p,u,y,z)  =  ^^  exp  j  -iky  {y+  z)  -[uiz-p)-  v(z-  x)]  /  c 


I  A  A  ^  A  A 

X  exp  iky  {y+  z)  ■[{p-x)  +  -{u-  v){y+  z)-{p  +  x)]  /  c 


(3.34) 


[  r  t+—T 

J  n  2 


^  1  ^ 

t - T 

V  2  y 


.kyiyyzHu-.),^^ 


J{p,u,y,z)d”'yd"z 


The  integral  in  square  brackets  is  the  radar  ambiguity  function 

^  '  1  V  r  1  ^ 


z(^,^)=  j 


=  I  ^  + 

V  / 


t  —  r 

V  2  y 


e"^‘dt 
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Hence,  the  imaging  PSF  can  be  written  as 


(  AAA  A 

i:( j!7, M , z)  =  JJ exp  -iky {y+z)-[u{y p)-v{z-x)\l c 


I  A  A  ^  A  A 

X  exp \  iky {y+  z)  ■[(p-x)  +  -(u- v)(y+  z)-(p  +  x)]  / c 


(3.35) 


X  J 


ky  {y+  z)-{u-  v),  {y+  z)-{p-x)/c  J{p,u,  z)d"'yd"z 


The  point  spread  i\mciion,K{p,u,y,z) ,  depends  on  the  geometry  of  the 
transmitter  and  receivers  and  also  on  the  transmitting  waveform  used.  When  it  is  shift 
invariant,  it  can  be  used  to  determine  the  position  and  velocity  resolution  of  the  imaging 
scene  for  a  given  imaging  configuration.  By  examining  K ,  one  can  effectively  model  the 
image  scene  with  meaningful  resolution. 
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IV.  RESULTS  AND  CONCLUSION 


This  chapter  describes  how  the  imaging  model  developed  in  Chapter  III  is 
implemented  and  analyzed.  Using  the  same  notion  described  in  Chapter  III,  three  sets  of 
MATLAB  programs  were  developed  for  the  analysis.  In  this  thesis,  a  four-parameter 
imaging  model  is  considered,  using  the  two-dimensions  of  position  jc  =  and  two- 

dimensions  of  velocity  v  =  (Vj,V2)  . 

A,  MODELS 

1.  The  Scattering  Model 

The  scattering  model  \j/ ,z ,t) ,  as  described  in  Equation  (3.19),  is  basically  a 
representation  of  the  measured  data  from  any  scatterers.  It  contains  information  of  the 
position  jc  =  (Xj,X2)  and  velocity  v  =  (Vj,V2)  of  the  unknown  scatterers.  A  MATLAB 

code  was  built  based  on  a  four-parameter  scattering  model  and  different  types  of 
waveforms. 

2,  The  Resolution  Model 

The  point  spread  function,  K{p,u,y,z)  ,  as  described  in  Equation  (3.35),  basically 
defines  the  resolution  of  the  imaging  scene.  It  depends  on  the  geometry  of  the 
transmitters  and  receivers  location  and  also  on  the  types  of  waveforms  used. 
Modifications  were  made  to  the  codes  developed  by  previous  theses  [10,11]  to  cater  to  a 
larger  number  of  transmitters  and  receivers. 
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3. 


The  Imaging  Model 


The  imaging  model,  as  described  in  Equation  (3.28),  is  the  main  focus  of  this 
thesis.  The  imaging  model  uses  the  cross  correlation  between  the  expected  target  position 
and  velocity  with  the  scattering  model  to  localize  the  position  and  velocity  of  the 
scatterer. 

The  focus  of  this  thesis  is  the  imaging  of  point  targets.  The  imaging  model  in 
Equation  (3.28)  is  not  restricted  to  any  type  of  transmitted  waveform  and  is  independent 
of  the  geometry  of  the  transmitter/receiver  configuration.  The  objective  of  this  thesis  is  to 
study  the  behavior  of  this  new  imaging  model  with  single  pulse  waveforms  in  different 
transmitter/receiver  configurations. 

Two  configurations  are  considered.  The  transmitter/receivers  geometries  used  for 
the  simulation  are  shown  in  Figure  4(a)  and  4(b).  The  receivers  are  denoted  by  “A”, 
while  the  transmitter  is  denoted  by  “x”.  Figure  4(a)  shows  the  geometry  of  the  linear 
array  of  receivers,  while  Figure  4(b)  shows  the  geometry  of  the  circular  array  of 
receivers. 


Linear  Array  Configuration 


Circular  Array  Configuration 


-10  -5  0  5  10 

X-Component  (km) 

(b) 


Figure  4.  Geometry  of  The  Transmitter  and  Receivers  in  (a)  Finear  Array 
Configuration  and  (b)  Circular  Array  Configuration 
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B. 


THE  SCATTERING  MODEL 


Assuming  the  simplest  ease,  in  whieh  there  is  one  transmitter  and  one  reeeiver, 
both  loeated  at  same  position,  10  km  away  from  the  origin  (-10km,  0),  we  eonsider  three 
targets  loeated  at  (-1000,  0),  (-200,0),  and  (400,0)  respeetively.  The  seattering  data  from 
the  three  targets  are  shown  in  Figure  5. 


Scanertng  Data  Of  Targets  At  Various  Locations 


Figure  5.  Seattering  Data  from  Targets  At  Various  Loeations 


From  the  scattering  data  shown  in  Figure  5.  ,  it  is  observed  that  the  scattering 
signal  is  a  time-delayed  version  of  the  transmitted  pulse.  The  further  away  the  target  is 
from  the  transmitter  and  receiver,  the  greater  the  time  delay. 

Figure  6.  shows  the  scattering  data  from  these  targets  using  the  linear  array 
configuration,  and  Figure  6.  shows  the  scattering  data  from  the  targets  using  the  circular 
array  configuration. 
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Linear  Array  Configuration 


Circular  Array  Configuration 


Figure  6.  Scattering  Data  From  Targets  At  Various  Locations 


C.  THE  RESOLUTION  MODEL 

The  position  and  velocity  resolutions  of  the  image  are  determined  by  their  point 
spread  function,  which  depends  on  the  pulse  width  and  the  transmitter  and  receiver 
geometries.  In  order  to  analyze  how  various  waveforms  and  geometries  affect  the 
resolution  of  the  imaging  system,  two  simple  waveforms  were  selected;  namely  a  single 
rectangular  pulse  and  a  single  chirp  pulse.  The  geometries  used  were  the  linear  array  and 
circular  array  configuration  shown  in  Figure  4.  . 

1,  Single  Rectangular  Pulse 

Using  a  single  rectangular  pulse  of  width  of  1  micro  second  as  the  transmitted 
waveform,  the  -3  dB  position  resolution  obtained  for  the  linear  array  configuration  and 
circular  array  configuration  is  shown  in  Figure  7.  .  The  -3  dB  velocity  resolution 

obtained  for  the  linear  array  configuration  and  circular  array  configuration  is  shown  in 
Figure  8.  . 
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Linear  Array  Configuration 


Circular  Array  Configuration 


Figure  7.  The  -3  dB  position  resolution  of  1  jus  reetangular  pulse  for  (a)  linear 
array  configuration  and  (b)  circular  array  configuration 


Linear  Array  Configuration 


-200  -too  0  100  200 

X-Component  (m) 


(a) 


Circular  Array  Configuration 


X-Component  (m) 


(b) 


Figure  8.  The  -3  dB  velocity  resolution  of  1  jus  rectangular  pulse  for  (a)  linear 
array  configuration  and  (b)  circular  array  configuration 


The  linear  array  configuration  has  a  position  resolution  that  is  narrower  in  the  x- 
components  and  broader  in  the  y-components,  while  the  position  resolution  has  a  more 
uniform  resolution  in  both  the  x-components  and  y  components.  However,  both 
configurations  have  very  poor  velocity  resolution  (which  is  not  meaningful  for  spatial 
imaging  purposes). 
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By  using  a  longer  rectangular  pulse  with  pulse-width  of  50  micro  seconds  as  the 
transmitted  waveform,  the  -3  dB  position  resolution  obtained  for  the  linear  array 
configuration  and  circular  array  configuration  is  shown  in  Figure  9.  .  The  -3  dB  velocity 
resolution  obtained  for  the  linear  array  configuration  and  circular  array  configuration  is 
shown  in  Figure  8.  . 


Linear  Array  Configuration  Circular  Array  Configuration 


X-Component  (m)  X-Component  (m) 


(a)  (b) 

Figure  9.  The  -3  dB  position  resolution  of  50  //s'  rectangular  pulse  for  (a)  linear 
array  configuration  and  (b)  circular  array  configuration 


Linear  Array  Configuration  Circuair  Array  Configuration 


Figure  10.  The  -3  dB  velocity  resolution  of  50  jus  rectangular  pulse  for  (a)  linear 

array  configuration  and  (b)  circular  array  configuration 

As  shown  in  Figure  9.  ,  both  configurations  have  poor  position  resolution. 

However,  as  shown  in  Figure  10.  the  velocity  resolution  has  improved. 
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These  examples  show  that  a  smaller  pulse-width  has  good  position  resolution, 
while  a  broader  pulse-width  has  good  velocity  resolution.  The  geometry  does  affect  the  x 
and  y  components  of  the  both  the  velocity  and  position  resolution. 


Single  Chirp  Pulse 


A  single  chirp  pulse  with  a  one  micro  second  pulse-width  having  a  chirp  rate  of 
0.5x10'^  Hz/s  is  used  as  the  transmitted  waveform.  The  -3dB  position  resolution 
obtained  for  the  linear  array  and  circular  array  is  shown  in  Figures  1 1  The  -3  dB  velocity 
resolution  obtained  for  the  linear  array  and  circular  array  is  shown  in  Figure  12. 


Linear  Array  Configuration 


Circular  Array  Configuration 


100  50 

X-Component  (m) 

(a) 


-200 


-150  -100  -50 

X-Component  (m) 

(b) 


Figure  1 1 .  The  -3  dB  position  resolution  of  1  /us  pulse  with  chirp  rate  of 

0.5x10'^  Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 
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X-Component  (m) 


X-Component  (m) 


(a) 


(b) 


Figure  12.  The  -3  dB  velocity  resolution  of  1  jus  pulse  with  chirp  rate  of 

0.5  X lO'^  Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 


From  Figures  11  and  12,  it  is  observed  that  the  position  and  the  velocity 
resolution  is  the  same  as  the  rectangular  pulse  of  1  micro  second  width.  Where  there  is 
good  position  resolution  and  poor  velocity  resolution. 

Figures  13  and  14  show  the  position  and  velocity  resolution  of  a  higher  chirp  rate. 
Using  a  higher  chirp  rate  of  lOxlO'^Hz/s,  the  position  resolution  has  improved 
significantly.  However,  the  velocity  resolution  is  still  poor. 


Linear  Array  Configuration  Circular  Array  Configuration 


Figure  13.  The  -3  dB  position  resolution  of  1  jus  pulse  with  chirp  rate  of 

10x10'^  Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 
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200 


Linear  Array  Configuration 


Circular  Array  Configuration 


(a) 


(b) 


Figure  14.  The  -3  dB  velocity  resolution  of  1  jus  pulse  with  chirp  rate  of 

lOxlO'^Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 


A  broader  pulse  width  is  used  to  analyze  the  effects  of  both  the  position  and 
velocity  resolution.  Using  a  fifty  micro  second  pulse  width  and  a  chirp  rate  of 
0.5  X  lO'^  Hz/s,  the  results  are  shown  in  Figures  15  and  16. 


Linear  Array  Configuration 


Circular  Array  Configuration 
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200 


Figure  15.  The  -3  dB  position  resolution  of  50  jis  pulse  with  chirp  rate  of 

0.5  X  lO'^  Hz/s  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 
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Linear  Array  Configuration 


Circular  Array  Configuration 


X-Component  (m)  X-Component  (m) 


(a)  (b) 

Figure  16.  The  -3  dB  velocity  resolution  of  50  /us  pulse  with  chirp  rate  of 

0.5x10'^  Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 


From  Figures  15  and  16,  it  is  observed  that  the  position  and  velocity  resolutions 
have  improved.  Figures  17  and  18  shows  the  results  of  a  5Q/us  pulse  width  and  a  chirp 

rate  of  lOxlO'^Hz/s. 


Linear  Array  Configuration  Circular  Array  Configuration 


X-Component  (m)  X-Component  (m) 


(a)  (b) 

Figure  17.  The  -3  dB  position  resolution  of  50  /us  pulse  with  chirp  rate  of 

10  X  lO'^  Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 
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Linear  Array  Configuration 


Circular  Array  Configuration 


(a)  (b) 

Figure  18.  The  -3  dB  velocity  resolution  of  50  /us  pulse  with  chirp  rate  of 

lOxlO'^Hz/s  for  (a)  linear  array  configuration  and  (b)  circular  array 

configuration 


From  Figure  17,  it  is  clearly  seen  that  the  position  resolution  is  very  fine,  which  is 
good  for  localizing  targets,  however  Figure  18  shows  no  improvement  in  the  velocity 
resolution. 

From  the  chirp  pulse  analysis,  it  is  observed  that  velocity  resolution  depends  on 
the  length  of  the  pulse-width;  the  broader  the  pulse,  the  finer  the  velocity  resolution.  The 
position  resolution  depends  on  the  chirp  rate,  the  higher  the  chirp  rate,  the  finer  the 
position  resolution.  Waveform  selection  is  important  for  imaging  moving  targets;  it  must 
have  the  ability  to  resolve  the  target’s  position  and  velocity  of  interest. 

D.  THE  IMAGING  MODEL 

Three  targets  are  used  to  analyze  the  behavior  of  the  imaging  model.  The  targets 
position  X  and  velocity  v  are  Xj(-35,25),  X2(-20,-30),  X3(-40,45),  and  Vj(0,-5), 

V2(10,10),  V3(-5,-10)  The  transmitter  and  receiver  geometries  are  the  linear  array 
configuration  and  circular  array  configuration  as  described  earlier. 
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Linear  Array  Configuration 


Figures  19,  20,  and  21  show  the  image  output  based  on  the  expected  velocity 
Mj(0,-5),  M2(10,10)  and  M3(-5,-10)  respectively. 


Image  Output  of  Expected  Velocity  u(0,-5)  Image  Output  of  Expected  Velocity  u(0.-5) 

100  r - ^ ^ H 


Figure  19.  Image  output  of  expected  velocity  of  Wj(0,-5) 


Image  Output  of  Expected  Velocity  u(10,10)  Image  Output  of  Expected  Velocity  u(10,10) 

100  r - ^ ^ ^ - n 


Figure  20.  Image  output  of  expected  velocity  of  W2(10,10) 
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Image  Output  of  Expected  Velocity  u(-5,-10)  Image  Output  of  Expected  Velocity  u(-5.-10) 
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Figure  21.  Image  output  of  expected  velocity  of  U3(-5,-10) 


All  the  images  output  show  that  the  target  positions  are  localized  accurately 
without  significant  artifacts.  A  threshold  level  of  3dB  lower  than  the  peak  is  used  to  reject 
the  image  artifacts.  Figure  22  shows  the  image  output  after  the  threshold  is  applied 


Image  Output  Image  output 

100  r - ^ ^ ^ 1 


Figure  22.  Threshold  image  output  of  all  expected  velocity  space 


Figure  22.  shows  a  very  clean  image  output  with  corresponding  velocity  vectors 
assigned.  This  is  achievable  mainly  due  to  the  cross  correlation  value  of  the  artifacts, 
which  are  small  compared  to  the  cross  correlation  values  that  are  “matched”. 
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Circular  Array  Configuration 


Figures  23,  24,  and  25  show  the  image  output  based  on  the  expected  velocity, 
Mj(0,-5),  M2(10,10)  and  M3  (-5, -10)  respectively 


Image  Output  of  Expected  Velocity  u(0,-5) 


Image  Output  of  Expected  Velocity  u(0.-5) 
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Figure  23.  Image  output  of  expected  velocity  of  Mj(0,-5) 


Image  Output  of  Expected  Velocity  u(10.10) 


1000 


m  500 


100 


X-Component  (m) 


Y-component  (m)  ’’'OO  -100  X-component  (m) 


Figure  24.  Image  output  of  expected  velocity  of  M2(10,10) 
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Figure  25. 


Image  output  of  expected  velocity  of  M3(-5,-10) 


All  the  images  show  that  the  target  positions  are  localized  accurately;  however,  in 
this  configuration  there  are  more  image  artifacts.  A  threshold  level  of  3dB  lower  than  the 
peak  is  used  to  reject  the  image  artifacts.  Figure  26  shows  the  image  output  after  the 
threshold  is  applied 


Image  Output  Image  Output 


Figure  26.  Image  output  of  all  expected  velocity  space 


The  three  targets  are  accurately  localized  in  both  position  and  velocity.  However, 
due  to  the  image  artifacts  as  shown  in  Figures  23,  24,  and  25,  the  -3  dB  threshold  level  is 
not  sufficient  to  reject  all  image  artifacts 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  analyzes  the  new  imaging  approaeh  developed  by  Cheney  and  Borden 
[1]  that  ean  aeeommodate  target  motion  during  the  imaging  proeess.  The  simulation 
results,  obtained  by  using  MATLAB,  showed  that  the  new  imaging  seheme  is  well 
behaved.  It  is  able  to  loealize  the  target  in  position  and  veloeity.  In  addition,  the  geometry 
of  the  transmitters,  reeeivers,  and  waveforms  utilized  affeet  the  behavior  of  the  imaging 
system. 

This  thesis  used  single  pulse  radar  waveforms,  but  in  reality,  fielded  radar  systems 
typieally  utilize  long  eoded  pulse  trains.  In  addition,  multistatie  radar  systems  may  not  be 
restrieted  to  one  transmitter.  In  this  ease,  the  imaging  seheme  has  to  be  adjusted  to 
ineorporate  more  realistie  radar  waveforms  and  optimized  for  the  transmitter/reeeiver 
geometries  for  implementation  of  the  imaging  algorithm.  These  will  be  used  to  support 
the  development  of  the  eventual  imaging  algorithm.  Finally,  real  world  target  data  ean 
then  be  applied  to  the  developed  imaging  algorithm  to  assess  its  performanee. 


42 


APPENDIX:  MATLAB  CODES 


Scattering  Data  of  Linear  Array  Configuration 

close  all 
clear  all 
clc 

%Origin  of  the  image  scene  is  set  at  (0,0) 

%Target  is  set  at  (1km,  1km)  away  from  the  origin  with  Velocity  of  (1,0) 

%Transmitted  is  set  at  (-10km,0)  away  from  the  origin 

%Target  Information 

tt=[-1000  0;-200  0;400  0];  %;0  840;50  50];  %Target  X  position,  Y  position 

tt_vel=[0  0;0  0;0  0];  %;0  170;800  -850];  %Target  velocity  in  X  direction,  Y  direction 
N_tt=length(tt(: ,  1 ));  %No  of  T argets 

%Transmitter  Information 

Tx  =[-I0e3  0];  %Transmitter  X  position,  Y  position 
N_Tx=length(Tx(:,I));  %No  of  Transmitter 

T_tx=0;  %Start  time  of  transmitted  pulse  =  always  zero  for  the  case  of  single 

transmitter 

Tx_mag=sqrt(Tx(  1 , 1  )^2-l-T  x(  1 ,2)^2); 

Tx_hat=Tx^qrt(Tx(  1 , 1  )^2-l-Tx(  1 ,2)^2); 

%Receiver  Information 

%single  line  of  51  receivers  at  X=10e3  away  from  the  origin. 

Rx_x(  1 :5 1 , 1  )=- 1 0e3 ; 

Rx_y=fliplr(-25000;1000:25000)’; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx(:,l).*Rx(:,l)+Rx(:,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 

%Signal  information 
c=3e8; 

W_tx=2*pi*(10e9);  %Carrier  Freq  is  lOGHz 
K_tx=W_tx/c;  %Ky=Wy/c 

%Waveform  information 
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%pulse  width  of  0.2  us  (30m  of  resolution  and  listening  time  of  lOOus  (15km  of  Runamb) 
%  total  period  lOOus 

fs  =  20e6;  %  Using  baseband  signal  of  IkHz.  Therefore  20MHz 

sampling.  Nyquisting  sampling  rate 
ts  =  1/fs; 

tl  =0:1 :20;  %0:ts:  le-6;  %0: 1 :20;  %pulse  transmit  time 

t2  =  0:1:1979;  %le-6+ts:ts:100e-6;  %0: 1:1979;  %listening  time 
period  =  lOOe-6; 

T_period=0:ts:  1  OOe-6; 

T=[tl  t2];  %period  of  lOOus 

s  =  [reotpuls(0,tl)  zeros(l,length(t2))];  %  0.2  us  Rect.  pulse 

S=fft(s); 

w=(2*pi/period)*T; 

TT_Data=zeros(  1,2001); 

%Generating  Target  Signal 
for  1=  1  :N_T x  %For  all  T ransmitter 

for  m=  1  :N_Rx  %F or  all  Receiver 

for  n=  1  :N_tt  %For  all  targets 

tau=  T_tx+(((Tx_mag(l,:)-(Tx_hat(l,:)*tt(n,:)'))+(Rx_mag(m,:)- 

(Rx_hat(m, :)  * tt(n, : )')))/c) ; %time  delay 
phi=K_tx*Rx_mag(m,:)- 

K_tx*(Tx_hat(l,:)+Rx_hat(m,:))*(tt(n,:)+((Rx_hat(m,:)*(Rx(m,:)-tt(n,:))')*tt_vel(n,:)/c))'; 

alpha=l-(Tx_hat(l,:)+Rx_hat(m,:))*(tt_vel(n,:)/c)'; 

TT_Data  =  TTData  +  exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 
end 

figure 

plot(T_period,abs(TT_Data)); 

Scattering  Data  of  Circular  Array  Configuration 

close  all 
clear  all 
clc 

%Origin  of  the  image  scene  is  set  at  (0,0) 

%Target  is  set  at  (1km,  1km)  away  from  the  origin  with  Velocity  of  (1,0) 

%Transmitted  is  set  at  (-10km,0)  away  from  the  origin 
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%Target  Information 

tt=[-1000  0;-200  0;400  0];  %;0  840;50  50];  %Target  X  position,  Y  position 

tt_vel=[0  0;0  0;0  0];  %;0  170;800  -850];  %Target  velocity  in  X  direction,  Y  direction 
N_tt=length(tt(: ,  1 ));  %No  of  T argets 

%Transmitter  Information 

Tx  =[-10e3  0];  %Transmitter  X  position,  Y  position 
N_Tx=length(Tx(:,l));  %No  of  Transmitter 

T_tx=0;  %Start  time  of  transmitted  pulse  =  always  zero  for  the  case  of  single 

transmitter 

T  x_mag=sqrt(Tx(  1 , 1  )^2-l-T  x(  1 ,2)^2); 

T  x_hat=T  x/ sqrt(Tx(  1 , 1  )^2-l-Tx(  1 ,2)^2); 

%Receiver  Information 

%Receiver  in  Circular  Configuration  at  10km  away  from  the  origin. 
theta=0:.125;2*pi; 

[Rx_x  Rx_y]=pol2cart(theta,10e3);  %R  is  1 0km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx(: ,  1 ).  *Rx( : ,  1  )+Rx(:  ,2).  *Rx( :  ,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

%Signal  information 
c=3e8; 

W_tx=2*pi*(10e9);  %Carrier  Freq  is  lOGHz 
K_tx=W_tx/c;  %Ky=Wy/c 

%Waveform  information 

%pulse  width  of  0.2  us  (30m  of  resolution  and  listening  time  of  lOOus  (15km  of  Runamb) 
%  total  period  lOOus 

fs  =  20e6;  %  Using  baseband  signal  of  IkHz.  Therefore  20MHz 

sampling.  Nyquisting  sampling  rate 
ts  =  1/fs; 

tl  =0:1 :20;  %0:ts:  le-6;  %0: 1 :20;  %pulse  transmit  time 

t2  =  0:1:1979;  %le-6+ts:ts:100e-6;  %0: 1:1979;  %listening  time 
period  =  lOOe-6; 

T_period=0:ts:  1  OOe-6; 

T=[tl  t2];  %period  of  lOOus 

s  =  [rectpuls(0,tl)  zeros(l,length(t2))];  %  0.2  us  Rect.  pulse 

S=fft(s); 

w=(2*pi/period)*T; 
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TT_Data=zeros(  1,2001); 

%Generating  Target  Signal 
for  1=  1  :N_T x  %For  all  T ransmitter 

for  m=l  ;N_Rx  %For  all  Receiver 

for  n=  1  :N_tt  %For  all  targets 

tau=  T_tx+(((Tx_mag(l,:)-(Tx_hat(l,:)*tt(n,:)'))+(Rx_mag(m,:)- 

(Rx_hat(m, :)  * tt(n, : )')))/c) ; %time  delay 
phi=K_tx*Rx_mag(m,:)- 

K_tx*(Tx_hat(l,:)+Rx_hat(m,:))*(tt(n,:)+((Rx_hat(m,:)*(Rx(m,:)-tt(n,:))')*tt_vel(n,:)/c))'; 
alpha=  1  -(T  x_hat(l,  :)+Rx_hat(m, :))  *  (tt_vel(n, :)/ c)'; 

TT_Data  =  TTData  +  exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 
end 

plot(T_period,abs(TT_Data)); 

Position  PSF  of  Linear  Array  Configuration  Using  Rectangular  Pulse 

clear  all 
close  all 
clc 

%%  Variables  to  be  specified 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%Transmitter 
Tx=[-10e3  0]; 

N_Tx=length(Tx(:,l)); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=Tx^qrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver 

Rx_x(  1  ;5 1 , 1  )=- 1 0e3 ; 

Rx_y=fliplr(-25000;1000;25000)'; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 
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E_tt_vel=[0  0]; 

%%  Parameters  Defined 
P_Width=  le-6; 

J=l;  %  Jaeobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%Expected  Target  Position 
res=l; 

E_tt_x=-50:res:50;  %Sampled  by  Range  Resolution 
E_tt_y=- 1  OOires:  1 00;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

value=0; 

for  g=l  ;length(E_tt_y) 
for  h=l  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

for  1=1  :N_Tx  %Por  all  Transmitter 
for  m=l  ;N_Rx  %Por  all  Receiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 
(tt_vel  *  (Rx_hat(m, : )  *  tf  )))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  ;)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_liat(l,  :)+Rx_hat(m, :))  *(E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=pi*(P_Width-abs(tau))*exp(-i*0.5*neu*P_Width)*sinc((0.5/pi)*(- 
neu)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 
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KT=K-0.707; 

K_3db=(KT>0).*K; 

figure 

imagesc(E_tt_x,E_tt_y,abs(k)) 

title('Position  PSE  of  Einear  Array  of  Receivers'); 

figure 

imagesc(E_tt_x,E_tt_y,K) 

colormap(gray); 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_x,E_tt_y  ,K_3  db) 
colormap(gray); 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

Velocity  PSF  of  Linear  Array  Configuration  Using  Rectangular  Pulse 

close  all 
clear  all 
clc 

%%  Variables  to  be  specified 

%  Actual  Target  Position  and  Velocity 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%  Transmitter  location 
Tx=[-10e3  0]; 

N_Tx=length(Tx(:,l)); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=Tx/sqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  Eocation 
Rx_x(  1  ;5 1 , 1  )=- 1 0e3 ; 

Rx_y=fliplr(-25000:1000:25000)'; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 
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E_tt=[0  0]; 

P_Width=le-6; 

%%  Parameters  Defined 
J=l;  %  Jaeobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%%  PSF  Velocity  Defined 
res=l; 

E_tt_vel_x=- 1 50:res:  150;  %Sampled  by  Range  Resolution 
E_tt_vel_y=- 1 50:res:  150;  %Sampled  by  Range  Resolution 
E_tt_vel_y=fliplr(E_tt_vel_y); 

value=0; 

for  g=l  :length(E_tt_vel_y) 
for  h=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,h)  E_tt_vel_y(  1  ,g)] ; 
for  1=1  :N_Tx  %For  all  Transmitter 
for  m=l  ;N_Rx  %For  all  Receiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 
(tt_vel  *  (Rx_hat(m, : )  *  tf  )))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  ;)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_liat(l,  :)+Rx_hat(m, :))  *(E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=pi*(P_Width-abs(tau))*exp(-i*0.5*neu*P_Width)*smc((0.5/pi)*(- 
neu)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 
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KT=K-0.99995; 

K_3db=(KT>0).*K; 

figure 

imagesc(E_tt_vel_x,E_tt_vel_y,abs(k)) 

titleC Velocity  PSE  of  Einear  Array  of  Receivers'); 

figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K) 

titleC Velocity  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K_3db) 

colormap(gray); 

titleC  Velocity  PSE  of  Einear  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

Position  PSF  of  Circular  Array  Configuration  Using  Rectangular  Pulse 

close  all 
clear  all 
clc 

%%  Variables  to  be  specified 

tt=[0  0];  %  actual  target 

tt_vel=[0  0];  %  actual  target  velocity 

%Transmitter 
Tx=[-10e3  0]; 

N_Tx=length(Tx(:,l)); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  in  Circular  Configuration  at  10km  away  from  the  origin. 
theta=0;.125:2*pi; 

[Rx_x  Rx_y]=pol2cart(theta,10e3);  %R  is  10km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 


%%  Parameters  Defined 
P_Width=le-6; 

J=l;  %  Jacobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%Expected  Target  Position 
res=l; 

E_tt_x=-200;res:200;  %Sampled  by  Range  Resolution 
E_tt_y=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

value=0; 

for  g=l  :length(E_tt_y) 
for  h=  1  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

for  1=1  :N_Tx  %Eor  all  Transmitter 
for  m=l  :N_Rx  %Eor  all  Receiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 
(tt_vel  *  (Rx_hat(m, : )  *  tf  )))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  ;)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=(P_Width-abs(tau))*exp(i*0.5*neu*(P_Width- 
abs(tau)))  *  sinc(0 . 5  *neu*  (P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 
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figure 

imagesc(E_tt_x,E_tt_y,abs(k)) 

title('Position  PSE  of  Circular  Array  of  Receivers'); 

figure 

imagesc(E_tt_x,E_tt_y,K) 

title('Position  PSE  of  Circular  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_x,E_tt_y  ,K_3  db) 
colormap(gray); 

title('Position  PSE  of  Circular  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Position  PSE  of  Circular  Array  of  Receivers  -  Normalised  @  3dB  Threshold 
Pulse  Width '  num2str(P_Width)]; 
saveas(gcf, filename, 'fig') 
close 

Velocity  PSF  of  Circular  Array  Configuration  Using  Rectangular  Pulse 


%  Actual  Target  Position  and  Velocity 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%  Transmitter  location 
Tx=[-10e3  0]; 

N_Tx=length(Tx(:,l)); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=Tx/sqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  in  Circular  Configuration  at  10km  away  from  the  origin. 
theta=0;.125:2*pi; 

[Rx_x  Rx_y]=pol2cart(theta,10e3);  %R  is  10km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx(: ,  1 ).  *Rx( : ,  1  )+Rx(:  ,2).  *Rx( :  ,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 
Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 

E_tt=[0  0]; 

%  P_Width=le-6; 

%  gamma=0.5el2; 
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%%  Parameters  Defined 
J=l;  %  Jaeobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%%  PSF  Velocity  Defined 
res=l; 

E_tt_vel_x=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_vel_y=-200;res:200;  %Sampled  by  Range  Resolution 
E_tt_vel_y=fliplr(E_tt_vel_y); 

value=0; 

for  g=l  :length(E_tt_vel_y) 
for  h=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,h)  E_tt_vel_y(  1  ,g)] ; 
for  1=1  ;N_Tx  %f  or  all  Transmitter 
for  m=l  ;N_Rx  %f  or  all  Receiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 
(tt_vel*  (Rx_hat(m, :)  *tf  )))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  :)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_liat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=(P_Width-abs(tau))*exp(i*0.5*(neu+gamma*tau)*(P_Width- 
abs(tau)))*smc(0.5*(neu+gamma*tau)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 
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figure 

imagesc(E_tt_vel_x,E_tt_vel_y,abs(k)) 
title('Position  PSE  of  Einear  Array  of  Receivers'); 

figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K) 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K_3db) 

colormap(gray); 

titleC Velocity  PSE  of  Circular  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Velocity  PSE  of  Circular  Array  of  Receivers  Pulse  Width'  num2str(P_Width) 
'-  Normalised  @  3dB  Threshold  Chirp  Rate'  num2str(gamma)]; 
saveas(gcf, filename, 'fig') 
close 

Position  PSF  of  Linear  Array  Configuration  Using  Chirp  Pulse 

%%  Variables  to  be  specified 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%Transmitter 
Tx=[-10e3  0]; 

N_T  x=length(T  x( : ,  1 )); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver 

Rx_x(  1 :5 1 , 1  )=- 1 0e3 ; 

Rx_y=fliplr(-25000:1000:25000)'; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( ; ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 

%%  Parameters  Defined 
E_tt_vel=[0  0]; 

%  P_Width=  le-6;  %Pulse  Width 
%  gamma  =  0.5el2;  %Chirp  Rate 
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J=l;  %  Jacobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%Expected  Target  Position 
res=l; 

E_tt_x=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_y=-200;res:200;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

value=0; 

for  g=l  :length(E_tt_y) 
for  h=l  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

for  1=1  :N_Tx  %Eor  all  Transmitter 
for  m=l  ;N_Rx  %Eor  all  Reeeiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 

(tt_vel*(Rx_hat(m,:)*tf)))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  :)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=(P_Width-abs(tau))*exp(i*0.5*(neu+gamma*tau)*(P_Width- 
abs(tau)))*sinc(0.5*(neu+gamma*tau)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 

figure 
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imagesc(E_tt_x,E_tt_y,abs(k)) 

title('Position  PSE  of  Einear  Array  of  Receivers'); 

figure 

imagesc(E_tt_x,E_tt_y,K) 

colormap(gray); 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_x,E_tt_y  ,K_3  db) 
colormap(gray); 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Position  PSE  of  Einear  Array  of  Receivers  Pulse  Width'  num2str(P_Width)  '- 
Normalised  @  3dB  Threshold  Chirp  Rate'  num2str(gamma)]; 
saveas(gcf, filename, 'fig') 
close 

Velocity  PSF  of  Linear  Array  Configuration  Using  Chirp  Pulse 

%  Actual  Target  Position  and  Velocity 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%  Transmitter  location 
Tx=[-10e3  0]; 

N_T  x=length(T  x( : ,  1 )); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  Eocation 
Rx_x(  1 :5 1 , 1  )=- 1 0e3 ; 

Rx_y=fliplr(-25000:1000:25000)'; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( ; ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 

%%  Parameters  Defined 
E_tt=[0  0]; 

%  P_Width=50e-6; 

%  gamma=10el2; 
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J=l;  %  Jacobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%%  PSF  Velocity  Defined 
res=l; 

E_tt_vel_x=-200:res:200;  %Sampled  by  Veloeity  Resolution 
E_tt_vel_y=-200:res:200;  %Sampled  by  Velocity  Resolution 
E_tt_vel_y=fliplr(E_tt_vel_y); 

value=0; 

for  g=l:length(E_tt_vel_y) 
for  h=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,h)  E_tt_vel_y(  1  ,g)] ; 
for  1=1  :N_Tx  %f  or  all  Transmitter 

for  m=l  ;N_Rx  %f  or  all  Reeeiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 

(tt_vel*(Rx_hat(m,:)*tf)))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  :)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=(P_Width-abs(tau))*exp(i*0.5*(neu+gamma*tau)*(P_Width- 
abs(tau)))*sinc(0.5*(neu+gamma*tau)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 

figure 

imagesc(E_tt_vel_x,E_tt_vel_y,abs(k)) 
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titleC Velocity  PSF  of  Linear  Array  of  Receivers'); 
figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K) 

titleC Velocity  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K_3db) 

colormap(gray); 

titleC  Velocity  PSE  of  Einear  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Velocity  PSE  of  Einear  Array  of  Receivers  Pulse  Width'  num2str(P_Width)  '- 
Normalised  @  3dB  Threshold  Chirp  Rate'  num2str(gamma)]; 
saveas(gcf, filename, 'fig') 
close 

Position  PSF  of  Circular  Array  Configuration  Using  Chirp  Pulse 

%%  Variables  to  be  specified 

tt=[0  0];  %  actual  target 

tt_vel=[0  0];  %  actual  target  velocity 

%Transmitter 
Tx=[-10e3  0]; 

N_Tx=length(Tx(:,l)); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  in  Circular  Configuration  at  10km  away  from  the  origin. 
theta=0:.125:2*pi; 

[Rx_x  Rx_y]=pol2cart(theta,10e3);  %R  is  1 0km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx(: ,  1 ).  *Rx( : ,  1  )+Rx(:  ,2).  *Rx( :  ,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

E_tt_vel=[0  0]; 

%%  Parameters  Defined 
%  P_Width=le-6; 

%  gamma=0.5el2; 
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J=l;  %  Jacobian 

c=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%Expected  Target  Position 
res=l; 

E_tt_x=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_y=-200;res:200;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

value=0; 

for  g=l  :length(E_tt_y) 
for  h=l  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

for  1=1  :N_Tx  %Eor  all  Transmitter 
for  m=l  ;N_Rx  %Eor  all  Reeeiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 

(tt_vel*(Rx_hat(m,:)*tf)))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  :)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y  *  ((Tx_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt_vel-tt_vel)') ; 
amb=(P_Width-abs(tau))*exp(i*0.5*(neu+gamma*tau)*(P_Width- 
abs(tau)))*sinc(0.5*(neu+gamma*tau)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 

figure 

imagesc(E_tt_x,E_tt_y,abs(k)) 
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title('Position  PSF  of  Circular  Array  of  Receivers'); 
figure 

imagese(E_tt_x,E_tt_y,K) 

title('Position  PSE  of  Cireular  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_x,E_tt_y  ,K_3  db) 
eolormap(gray); 

title('Position  PSE  of  Cireular  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Position  PSE  of  Circular  Array  of  Receivers  Pulse  Width'  num2str(P_Width) 
'-  Normalised  @  3dB  Threshold  Chirp  Rate'  num2str(gamma)]; 
saveas(gef, filename, 'fig') 
elose 

Velocity  PSF  of  Circular  Array  Configuration  Using  Chirp  Pulse 

%%  Variables  to  be  specified 

%  Aetual  Target  Position  and  Veloeity 

tt=[0  0];  %  actual  target  -  tt 

tt_vel=[0  0];  %  actual  target  velocity  -  tt_vel 

%  Transmitter  loeation 
Tx=[-10e3  0]; 

N_T  x=length(T  x( : ,  1 )); 

Tx_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  in  Circular  Configuration  at  10km  away  from  the  origin. 
theta=0:.125;2*pi; 

[Rx_x  Rx_y]=pol2cart(theta,10e3);  %R  is  10km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *Rx( :  ,2)) ; 

Rx_hat(:,2)=  Rx(:,2)./sqrt(Rx(:,l).*Rx(:,l)-l-Rx(:,2).*Rx(:,2)); 

E_tt=[0  0]; 

%  P_Width=le-6; 

%  gamma=0.5el2; 
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%%  Parameters  Defined 
J=l;  %  Jaeobian 

e=3*10^8; 

k_y=2*pi*(10e9)/c;  %  wave  vector  @f=10  GHz 

%%  PSF  Velocity  Defined 
res=l; 

E_tt_vel_x=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_vel_y=-200:res:200;  %Sampled  by  Range  Resolution 
E_tt_vel_y=fliplr(E_tt_vel_y); 

value=0; 

for  g=l  :length(E_tt_vel_y) 
for  h=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,h)  E_tt_vel_y(  1  ,g)] ; 
for  1=1  ;N_Tx  %f  or  all  Transmitter 
for  m=l  ;N_Rx  %f  or  all  Receiver 

argl=(Tx_hat(l,:)+Rx_hat(m,:))*((E_tt_vel*(Rx_hat(m,:)*E_tf))- 
(tt_vel  *  (Rx_hat(m, : )  *  tf  )))'; 

phasel=exp(-i*k_y*argl/c); 

arg2=(T  x_hat(l,  :)+Rx_hat(m,  :))*  ((E_tt-tt)+(0 . 5  *  (E_tt_vel- 
tt_vel)*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt+tt)')))'; 
phase2=exp(i*k_y*arg2/c); 
tau=((T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt-tt)')/ c ; 
if  abs(tau)<P_Width 

neu=  k_y*((Tx_hat(l,:)+Rx_hat(m,:))*(E_tt_vel-tt_vel)'); 
amb=(P_Width-abs(tau))*exp(i*0.5*(neu+gamma*tau)*(P_Width- 
abs(tau)))*sinc(0.5*(neu+gamma*tau)*(P_Width-abs(tau))); 
else 
amb=0; 
end 

value=  value  +  phasel*phase2*amb*J; 
end 
end 

k(g,h)=value.^2; 

value=0; 

end 

end 

K=abs(k)/max(max(abs(k))); 

KT=K-0.707; 

K_3db=(KT>0).*K; 

figure 
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imagesc(E_tt_vel_x,E_tt_vel_y,abs(k)) 
title('Position  PSE  of  Einear  Array  of  Receivers'); 

figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K) 

title('Position  PSE  of  Einear  Array  of  Receivers  -  Normalised'); 
figure 

imagesc(E_tt_vel_x,E_tt_vel_y,K_3db) 

colormap(gray); 

titleC Velocity  PSE  of  Circular  Array  of  Receivers  -  Normalised  @  3dB  Threshold'); 

lilename=['Velocity  PSE  of  Circular  Array  of  Receivers  Pulse  Width'  num2str(P_Width) 
'-  Normalised  @  3dB  Threshold  Chirp  Rate'  num2str(gamma)]; 
saveas(gcf, filename, 'fig') 
close 

Imaging:  Linear  Array  Configurations. 

close  all 
clear  all 
clc 

%Origin  of  the  image  scene  is  set  at  (0,0) 

%Target  is  set  at  (1km,  1km)  away  from  the  origin  with  Velocity  of  (1,0) 

%Transmitted  is  set  at  (-10km,0)  away  from  the  origin 

%Target  Information 

tt=[-35  25;-20  -30;40  45];  %;0  840;50  50];  %Target  X  position,  Y  position 

tt_vel=[0  -5;10  10;-5  -10];  %;0  170;800  -850];  %Target  velocity  in  X  direction,  Y 
direction 

N_tt=length(tt(: ,  1 ));  %No  of  T argets 

%Transmitter  Information 

Tx  =[-I0e3  0];  %Transmitter  X  position,  Y  position 
N_Tx=length(Tx(:,I));  %No  of  Transmitter 

T_tx=0;  %Start  time  of  transmitted  pulse  =  always  zero  for  the  case  of  single 

transmitter 

Tx_mag=sqrt(Tx(  1 , 1  )^2-l-T  x(  1 ,2)^2); 

Tx_hat=TxAqrt(Tx(  1 , 1  )^2-l-Tx(  1 ,2)^2); 

%Receiver  Information 

%single  line  of  51  receivers  at  X=10e3  away  from  the  origin. 

Rx_x(  1 :5 1 , 1  )=- 1 0e3 ; 
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Rx_y=fliplr(-25000:1000:25000)’; 

Rx=[Rx_x  Rx_y];  %Receiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx(: ,  1 ).  *Rx( : ,  1  )+Rx(:  ,2).  *Rx( :  ,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *  Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

%Signal  information 
c=3e8; 

W_tx=2*pi*(10e9);  %Carrier  Freq  is  lOGHz 
K_tx=W_tx/c;  %Ky=Wy/c 

%Waveform  information 

%pulse  width  of  0.2  us  (30m  of  resolution  and  listening  time  of  lOOus  (15km  of  Runamb) 
%  total  period  lOOus 

fs  =  20e6;  %  Using  baseband  signal  of  IkHz.  Therefore  20MHz 

sampling.  Nyquisting  sampling  rate 
ts  =  1/fs; 

tl  =0:1 :20;  %0:ts:  le-6;  %0;  1 :20;  %pulse  transmit  time 

t2  =  0:1:1979;  %le-6+ts:ts:100e-6;  %0: 1:1979;  %listening  time 
period  =  lOOe-6; 

T_period=0:ts:  1  OOe-6; 

T=[tl  t2];  %period  of  lOOus 

s  =  [reotpuls(0,tl)  zeros(l,length(t2))];  %  0.2  us  Rect.  pulse 

S=fft(s); 

w=(2*pi/period)*T; 

TT_Data=zeros(  1,2001); 

%Generating  Target  Signal 
for  1=  1  :N_T x  %For  all  T ransmitter 

for  m=  1  :N_Rx  %F or  all  Reeeiver 

for  n=  1  :N_tt  %For  all  targets 

tau=  T_tx-l-(((Tx_mag(l,:)-(Tx_hat(l,:)*tt(n,:)'))-l-(Rx_mag(m,:)- 

(Rx_hat(m, :)  * tt(n, : )')))/c) ; %time  delay 
phi=K_tx*Rx_mag(m,:)- 

K_tx*(Tx_hat(l,:)-l-Rx_hat(m,:))*(tt(n,:)+((Rx_hat(m,:)*(Rx(m,:)-tt(n,:))')*tt_vel(n,:)/c))'; 
alpha=  1  -(T  x_hat(l,  :)-l-Rx_hat(m, :))  *  (tt_vel(n, :)/ c)'; 

TT_Data  =  TTData  +  exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 
end 
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%  figure 

%  plot(T_period,abs(TT_Data)); 

%Expected  Target  Position 

E_tt_x=- 1 00:5 : 1 00;  %Sampled  by  Range  Resolution 
E_tt_y=- 1 00:5 ;  1 00;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

Expected  Target  Velocity 
E_tt_vel_x=-50: 1 :50; 

E_tt_vel_y=-50: 1 :50; 

E_tt_vel_y=fliplr(E_tt_vel_y); 

E=length(E_tt_y)  *  length(E_tt_x) ; 

%Generating  Information  Database 
[X,Y]=meshgrid(E_tt_x,E_tt_y); 

Vel_El=zeros(length(E_tt_x),length(E_tt_y)); 

Vel_V=zeros(length(E_tt_x),length(E_tt_y)); 

Vel_W=zeros(length(E_tt_x),length(E_tt_y)); 

nn=l; 

This  for  loop  is  for  spanning  across  vested  velocity  space 
for  u=l  :length(E_tt_vel_y) 
for  v=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,v)  E_tt_vel_y(  1  ,u)] ; 

E_tt_vel_mag=sqrt(E_tt_vel(  1 , 1  )^2+E_tt_vel(  1 ,2)^2); 

%Generating  Expected  Target  Database 
for  g=l:length(E_tt_y) 
for  h=l  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

E_TT_Data=zeros(  1  ,length(T)); 
for  1=  1  :N_T x  %f  or  all  T ransmitter 

for  m=l  :N_Rx  %f  or  all  Receiver 

tau=  T_tx+(((Tx_mag(l,:)-(Tx_hat(l,:)*E_tt(l,:)'))+(Rx_mag(m,:)- 

(Rx_hat(m,:)*E_tt(l,:)')))/c);  %time  delay 
phi=K_tx*Rx_mag(m,  :)- 

K_tx  *  (T  x_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt(  1 , :  )+((Rx_hat(m,  :)*(Rx(m,:)- 
E_tt(l,:))')*E_tt_vel/c))'; 

alpha=  1  -(T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt_vel/ c)'; 

ETTData  =  ETTData  + 

exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 
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I(g,h)=E_TT_Data*TT_Data'; 

end 

end 

I_abs=abs(I); 

I_norm=I_abs/max(max(I_abs)); 

I_tt=(I_nomi-0.707); 

I_tt=ceil(I_tt)*E_tt_vel_mag; 

[row,col]=find(I_tt); 

Vel_Er(row,col)=E_tt_vel(  1,1); 

Vel_V(row,col)=E_tt_vel(  1 ,2); 
ifnn==l 

l_tt_sum=l_tt; 

else 

l_tt_sum=l_tt_sum+l_tt; 

end 

nn=nn+l; 

end 

end 

figure 

surf(E_tt_x,E_tt_y,l_tt_sum);  %or  l_norm*E_tt_vel_mag 
%title(['veloeity  Vector  [',  num2str(E_tt_vel), ']']); 
colorbar; 
hold  on 

quiverS  (X,  Y,l_tt_sum,  Vel_U,  Vel_V,  Vel_W); 

Imaging:  Circular  Array  Configurations. 

close  all 
clear  all 
clc 

%Origin  of  the  image  scene  is  set  at  (0,0) 

%Target  is  set  at  (1km,  1km)  away  from  the  origin  with  Velocity  of  (1,0) 

%Transmitted  is  set  at  (-10km,0)  away  from  the  origin 

%Target  Information 

tt=[-35  25;-20  -30;40  45];  %;0  840;50  50];  %Target  X  position,  Y  position 

tt_vel=[0  -5;10  10;-5  -10];  %;0  170;800  -850];  %Target  velocity  in  X  direction,  Y 
direction 

N_tt=length(tt(: ,  1 ));  %No  of  T argets 

%Transmitter  Information 
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Tx  =[-10e3  0];  %Transmitter  X  position,  Y  position 
N_Tx=length(Tx(:,l));  %No  of  Transmitter 

T_tx=0;  %Start  time  of  transmitted  pulse  =  always  zero  for  the  ease  of  single 

transmitter 

T  x_mag=sqrt(Tx(  1 , 1  )^2+T  x(  1 ,2)^2); 

T  x_hat=T  x/ sqrt(Tx(  1 , 1  )^2+Tx(  1 ,2)^2); 

%Receiver  Information 

%Reoeiver  in  Cireular  Configuration  at  10km  away  from  the  origin. 
theta=0:.125;2*pi; 

[Rx_x  Rx_y]=pol2oart(theta,10e3);  %R  is  10km  away  from  origin 
Rx=[Rx_x'  Rx_y'];  %Reoeiver  position  in  X,Y  co-ordinates 
N_Rx=length(Rx(:,l)); 

Rx_mag=  sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

Rx_hat( : ,  1  )=Rx( : ,  1 )  ./sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )-l-Rx( :  ,2) .  *Rx( :  ,2)) ; 

Rx_hat(  :,2)=Rx(:,2)./ sqrt(Rx( : ,  1 ) .  *  Rx( : ,  1  )+Rx(  :,2).*Rx(:,2)); 

%Signal  information 
c=3e8; 

W_tx=2*pi*(10e9);  %Carrier  Freq  is  lOGHz 
K_tx=W_tx/c;  %Ky=Wy/c 

%  Waveform  information 

%pulse  width  of  0.2  us  (30m  of  resolution  and  listening  time  of  lOOus  (15km  of  Runamb) 
%  total  period  lOOus 

fs  =  20e6;  %  Using  baseband  signal  of  IkHz.  Therefore  20MHz 

sampling.  Nyquisting  sampling  rate 
ts  =  1/fs; 

tl  =0:1 :20;  %0:ts:  le-6;  %0: 1 :20;  %pulse  transmit  time 

t2  =  0:1:1979;  %le-6+ts:ts:100e-6;  %0: 1:1979;  %listening  time 
period  =  lOOe-6; 

T_period=0:ts:  1  OOe-6; 

T=[tl  t2];  %period  of  lOOus 

s  =  [rectpuls(0,tl)  zeros(l,length(t2))];  %  0.2  us  Rect.  pulse 

S=fft(s); 

w=(2*pi/period)*T; 

TT_Data=zeros(  1,2001); 

%Generating  Target  Signal 
for  1=  1  :N_T x  %For  all  T ransmitter 

form=l:N  Rx  %For  all  Receiver 
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for  n=  1  :N_tt  %For  all  targets 

tau=  T_tx+(((Tx_mag(l,:)-(Tx_hat(l,:)*tt(n,:)'))+(Rx_mag(m,:)- 

(Rx_hat(m, :)  * tt(n, : )')))/c) ; %time  delay 
phi=K_tx*Rx_mag(m,:)- 

K_tx*(Tx_hat(l,:)+Rx_hat(m,:))*(tt(n,:)+((Rx_hat(m,:)*(Rx(m,:)-tt(n,:))')*tt_vel(n,:)/c))'; 
alpha=  1 -(Tx_hat(l,  :)+Rx_hat(m, :))  *  (tt_vel(n,  :)/e)'; 

TT_Data  =  TTData  +  exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 
end 

%  figure 

%  plot(T_period,abs(TT_Data)); 

%Expected  Target  Position 

E_tt_x=- 1 00:5 ;  1 00;  %Sampled  by  Range  Resolution 
E_tt_y=- 1 00:5 : 1 00;  %Sampled  by  Range  Resolution 
E_tt_y=fliplr(E_tt_y); 

%Expected  Target  Veloeity 
E_tt_vel_x=-50: 1 :50; 

E_tt_vel_y=-50: 1 :50; 

E_tt_vel_y=fliplr(E_tt_vel_y); 

E=length(E_tt_y)  *  length(E_tt_x) ; 

%Generating  Information  Database 
[X,Y]=meshgrid(E_tt_x,E_tt_y); 

Vel_El=zeros(length(E_tt_x),length(E_tt_y)); 

Vel_V=zeros(length(E_tt_x),length(E_tt_y)); 

Vel_W=zeros(length(E_tt_x),length(E_tt_y)); 

nn=l; 

This  for  loop  is  for  spanning  across  vested  velocity  space 
for  u=  1  :length(E_tt_vel_y) 
for  v=l  :length(E_tt_vel_x) 

E_tt_vel=[E_tt_vel_x(  1  ,v)  E_tt_vel_y(  1  ,u)] ; 

E_tt_vel_mag=sqrt(E_tt_vel(  1 , 1  )^2+E_tt_vel(  1 ,2)^2); 

%Generating  Expected  Target  Database 
for  g=l  :length(E_tt_y) 
for  h=  1  :length(E_tt_x) 

E_tt=[E_tt_x(l,h)  E_tt_y(l,g)]; 

E_TT_Data=zeros(  1  ,length(T)); 
for  1=  1  :N_T x  %Eor  all  T ransmitter 

form=l:N  Rx  %Eor  all  Receiver 
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tau=  T_tx+(((Tx_mag(l,:)-(Tx_hat(l,:)*E_tt(l,:)'))+(Rx_mag(m,:)- 

(Rx_hat(m,:)*E_tt(l,:)')))/c);  %time  delay 
phi=K_tx*Rx_mag(m,  :)- 

K_tx  *  (T  x_hat(l, :  )+Rx_hat(m, : ))  *  (E_tt(  1 , :  )+((Rx_hat(m,  :)*(Rx(m,:)- 
E_tt(l,:))')*E_tt_vel/c))'; 

alpha=  1  -(T  x_hat(l,  :)+Rx_hat(m, :))  *  (E_tt_vel/ c)'; 

E_TT_Data  =  ETTData  + 

exp(i*phi)*exp(i*W_tx*alpha*T).*ifft(S.*exp(-i*w*tau)); 
end 
end 

I(g,h)=E_TT_Data*TT_Data'; 

end 

end 

I_abs=abs(I); 

I_norm=I_abs/max(max(I_abs)); 

I_tt=(I_nomi-0.707);  %applying  -3dB  thresholding 
I_tt=eeil(I_tt)*E_tt_vel_mag; 

[row,eol]=find(I_tt); 

Vel_El(row,eol)=E_tt_vel(  1,1); 

Vel_V(row,eol)=E_tt_vel(  1 ,2); 
ifnn==l 

I_tt_sum=I_tt; 

else 

I_tt_sum=I_tt_sum+I_tt; 

end 

nn=nn+l; 

end 

end 

figure 

surf(E_tt_x,E_tt_y,I_tt_sum);  %or  I_norm*E_tt_vel_mag 
%title(['veloeity  Veetor  [',  num2str(E_tt_vel), ']']); 
eolorbar; 
hold  on 

quiverS  (X,  Y,I_tt_sum,  Vel_El,  Vel_V,  Vel_W); 
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